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Abstract 

We examine the connection between the nonhnear integral equation (NLIE) de- 
rived from hght-cone lattice and sine-Gordon quantum field theory, considered as a 
perturbed c = 1 conformal field theory. After clarifying some delicate points of the 
NLIE deduction from the lattice, we compare both analytic and numerical predictions 
of the NLIE to previously known results in sine-Gordon theory. To provide the basis 
for the numerical comparison we use data from Truncated Conformal Space method. 
Together with results from analysis of infrared and ultraviolet asymptotics, we find 
evidence that it is necessary to change the rule of quantization proposed by Destri 
and de Vega to a new one which includes as a special case that of Fioravanti et al. 
This way we find strong evidence for the validity of the NLIE as a description of the 
finite size effects of sine-Gordon theory. 
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1 Introduction 



In the theory of 1 + 1 dimensional integrable models, the investigation of the spectrum on a 
spacetime cylinder with finite spatial volume plays a very important role in exploring finite 
size effects and provides a possibility to determine many important physical characteristics. 

One of the first approaches to compute the finite volume spectrum was the the Ther- 
modynamic Bethe Ansatz (TBA) which was used to calculate the vacuum (Casimir) 
energy f^. The method was later extended to include ground states of charged sectors [Q. 
More recently, using analytic properties of the TBA equations extended for complex values 
of the volume parameter, an approach to get excited states was proposed in 0, ||, The 
method of Q to get excited states sheds light on the analytic structure of the dependence 
of scaling functions on the spatial volume and is the only method developed so far to deal 
with excited states in perturbations of RCFTs. Its main drawbacks are that (1) it can be 
used only for systems for which a TBA equation describing the vacuum in finite volume is 
known; (2) to obtain the equation for a given excited state one has to do analytic continua- 
tion for each case separately, and a major part of this continuation can only be carried out 
numerically. Because of the requirement of the knowledge of the vacuum TBA equation 
and the complications of the analytic continuation, this method is limited at present to 
integrable perturbations of Virasoro minimal models and other simple perturbed conformal 
field theories. 

Another method developed to study the finite volume energy levels of 1 + 1 dimensional 
field theories is the Truncated Conformal Space (TCS) approach 0. Its big advantage is 
that it makes possible to study the finite volume spectrum of any theory for which a 
description in terms of perturbed conformal field theory is known (integrability is not 
required). In practice, so far, it has been mainly used to study perturbations of Virasoro 
minimal models. In this paper we use an extension of the TCS method to perturbations 
of c = 1 conformal field theories in order to get a basis for comparison with the results 
obtained from the Non-Linear Integral Equation (NLIE) method to be described below. 
The fundamental limitation of this approach is that it only allows one to get numerical 
predictions for the energy levels. 

The method which is in the focus of this paper is that of the NLIE deduced from 
integrable light cone lattice regularisation. The light cone lattice model was proposed 
by Destri and de Vega in ||^ where a fermionic operator was constructed on the lattice 
which was shown to satisfy the Thirring equations of motion in the continuum limit. 
Therefore it was plausible to conjecture that the light cone lattice has something to do 
with a regularisation of the sine-Gordon/massive Thirring (sG/mTh) field theory. The 
states of the lattice model can be described in a Bethe Ansatz framework which can be 
recast in the form of a nonlinear integral equation using residue tricks. The equation 
describing the vacuum has been derived in P, |lOl and it was shown that in the ultraviolet 
limit it reproduces the correct value of the central charge c = 1. We remark that similar 
methods were independently introduced in Condensed Matter Physics by other authors 
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The first extension to describe the excited state spectrum was introduced by Fiora- 
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vanti et al. [0] where the spectrum of states containing only sohtons (and no antisoh- 
tons/breathers) has been described. An extension to generic excited states of even topo- 
logical charge was described by Destri and de Vega in ||13[ which however gave a different 
prescription from that of Fioravanti et al. when specialized to the multi-soliton states 
examined in ||12|. The contradiction was resolved in ||14| where it was shown, using the 
ultraviolet conformal weights calculated from the NLIE and a numerical comparison to 
TCS results, that the prescription given by Destri and de Vega leads to states which are 
not contained in the Hilbert space while the one proposed by Fioravanti et al. describes 
the spectrum of pure soliton states correctly. 

The NLIE treated in this paper corresponds in algebraic terms to the nontwisted afiine 
Kac-Moody algebra A^^^ = sl{2). An extension of the approach to simply-laced algebras 



of ADE type was first given in [T^ for the case of the vacuum. More recently, in [|T^] there 
appeared the extension to describe the spectrum of excited states. 

In this paper we describe a generalization of our previous results in to generic 
excited states with even value of the topological charge, by extending the framework pro- 
posed by Destri and de Vega in In doing so we show that their form of the NLIE 
is not entirely correct and give a new derivation of the NLIE which leads to the correct 
result and avoids certain subtle points which are present in all of the derivations known 
up to now. We examine analytically both the infrared and ultraviolet asymptotics of the 
scaling functions and show that the quantisation prescription can be chosen in a way which 
is consistent with the ultraviolet limit of the sG/mTh theory. We also give comparisons of 
the scaling functions coming from the NLIE to the ones predicted by the TCS approach. 

The major motivation for studying the NLIE is that it gives a method of describing all 
excited states in a single framework for all values of the (sine-Gordon) coupling constant, 
in contrast to the TBA approach. We will return to a detailed comparison of the two 
methods in the conclusions. 

2 The sine- Gordon/massive Thirring model as a per- 
turbed CFT 

2.1 c=l CFT 

To fix our conventions and to define certain objects which are used later, we give a brief 
summary of the c = 1 free boson with a target space of a circle of radius R, which describes 
the UV limit of the theory. The Lagrangian of this CFT is taken to be 



C = — [ d.^{x, t)d''^{x, t)dx , X G [0, L] , (2.1) 

where L is the spatial volume (i.e. the theory is defined on a cylindrical spacetime with 
circumference L). In the sequel we will also use the complex Euclidean coordinates z = 
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g27r(t-ix)/L^ ^ ^ ^2^{t+ix)/L_ classify the superselection sectors by the C/(l)^ x C/(l)^ 
Kac-Moody symmetry algebra, generated by the currents 

J{z) — idz(fi , J{z) — idz(fi . 
The left/right moving energy-momentum tensor is given by 



T{z) = j^dz^d,^ = J2 ^kz-'-' , f{z) = ^d,^d,^ = LkZ 

k=—oo k=—oo 
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The coefficients L„ and L„ of the Laurent expansion of these fields generate two mutually 
commuting Virasoro algebras. If we require the (quasi)periodic boundary conditions 

ip{x + L, t) = (p{x, t) + 27rmi? , m e Z , 

then the sectors are labelled by a pair of integers (n, m), where ^ is the eigenvalue of the 
total field momentum ttq 



TTo = y 7r{x, t)dx , 7r(a;, t) = -^dt^p{x, t) , 



and m is the winding number, i.e. the eigenvalue of the topological charge Q defined by 



L 



Q = j d:cV{x, t)dx . 







In the sector with quantum numbers (n, m), the scalar field is expanded in modes as 
follows: 

(p{x,t)= (p{z) + (t){z) , 
Hz) = ^^o-ip+'^ogz + i^-akZ''' , 

fc^O 

4>(z) = 2*^° ~ log ^ + ^ X] k^''^''' ' 

k^Q 

where the left and right moving field momenta p± (which are in fact the two U{1) Kac- 
Moody charges) are given by 

P± = ;^±2^^- (2.2) 
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The Virasoro generators take the form 



^ oo ^ oo 

fc=— oo fc=— oo 

where the colons denote the usual normal ordering, according to which the oscillator with 
the larger index is put to the right. 

The ground states of the different sectors (n, m) are created from the vacuum by the 
(Kac-Moody) primary fields, which are vertex operators of the form 



V(n,m)iz,z) =: expi{p+(f){z) + p_(i){z)) : . (2.3) 

The left and right conformal weights of the field V(n,m) (i-e. the eigenvalues of Lq and Lq) 
are given by the formulae 

A± = f . (2.4) 

The Hilbert space of the theory is given by the direct sum of the Fock modules built over 
the states 



\n, m) = V(n,^)(0, 0) \vac) , (2.5) 
with the help of the creation operators a_fe , a_fc A; > 0: 

Ti^ ^ {a_fei . ..a^kpa-h ■ ■■a-ijn, m), h,... kp, k, . . . lq e Z+} 

(n, m) 

The Hamiltonian is expressed in terms of the Virasoro operators as 

HcFT^-r [Lo + Lo-—] , (2.6) 



L V " ■ " 12 

where the central charge is c = 1. The generator of spatial translations is given by 

27r 

P = — (Lo - Lo) . (2.7) 
The operator Lq — Lq is the conformal spin which has eigenvalue nm on the primary field 

V{n, m) ■ 
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One can also introduce twisted sectors using the operator T of the spatial translation 
with one period x x+L. The primary fields V(„,m) as defined above satisfy the periodicity 
condition TV(„ m) = V(„ m)- If we require the more general twisted boundary condition 
labelled by a real parameter u 

T^V{n^rn) = exp (iuQ) V(n,m) , 

then we can generate superselection sectors for which n G Z + The only twisted 
boundary condition occurring in this paper will be the one with i/ = tt, which is necessary 
to describe the fermions of the massive Thirring model. 

The classification of the (modular invariant) c = 1 CFTs can be found in ||T8| and is 
not need in the sequel. We stress only that a particular c = 1 CFT is specified by giving 
the spectrum of the quantum numbers (n, m) such that the corresponding set of vertex 
operators (and their descendants) forms a closed and local operator algebra. The locality 
requirement is equivalent to the fact that the operator product expansions of any two such 
local operators is single valued in the complex plane of z. This condition, which is weaker 
than the modular invariance of the CFT, is the adequate one since we consider the theory 
on a spacetime cylinder and do not wish to define it on higher genus surfaces. 

2.2 Sine-Gordon/massive Thirring field theory 

The Lagrangian of sine-Gordon theory is given by 

where $ denotes a real scalar field, while that of the massive Thirring theory is of the 
following form: 

CmTh = j [^{ilud" + mo)^ - 1^7^^^%.^) dx , (2.9) 

describing a current-current selfinteraction of a Dirac fermion It is known that the two 
theories are deeply related provided their coupling constants satisfy 

^ _ 1 
47r 1 + gj-K 

Both models can be considered as the perturbations of a c = 1 CFT by a potential V : 

HsG/mTh = HcFT + V , V = X [ {V^^,o){z,z) + V(^_^fl){z,z)) dx , (2.10) 
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which is related to the bosonic Lagrangian (^) by the following redefinitions of the field 
and the parameters: 



(^ = V47r$, R=^, A = ^. (2.11) 



For later convenience, we also define a new parameter p with 



In terms of p = 1 is the free fermion point, and p = y , k = 1,2, . . . are the thresholds 

k 

where a new bound state (sG breather) appears, p < 1 corresponds to the attractive and 
p > 1 to the repulsive regime. The potential term becomes marginal when = Svr which 
corresponds to p = oo. The perturbation conserves the topological charge Q, which can be 
identified with the usual topological charge of the sG theory and with the fermion number 
of the mTh model. 

The difference between the two theories is that they correspond to the perturbation of 
two dijferent local c=l CFTs (but by the same operator). The short distance behaviour of 
the sG theory is described by the c=l CFT with a local operator algebra At generated by 
the vertex operators 

{V(n,m) ■ n, meZ}, 

all of which satisfy the periodic boundary condition and lead to a modular invariant par- 
tition function. The primary fields of the UV limit of mTh theory are 

{V(n,m) ■ n eZ, m e2Zorn eZ+^, m e2Z+l} , 

generating an operator algebra Af, which also contains fields having antiperiodic boundary 
conditions (i/ = tt and odd value of the topological charge), and corresponds to a 
invariant partition function. In particular, the fields V(±i/2,±i) have conformal spin | and 
describe the UV limit of the one-fermion states. The interested reader can find a detailed 
exposition in |0. Note that the two algebras share a common subspace with even values of 



the topological charge, generated by {V(n,m) '■ n G Z, m G 2Z}, where the massive theories 



described by the Lagrangians (|2.8D and (|2^) are identical. It is exactly this subspace that 



is accessible with the NLIE technique and is discussed in this paper. 



3 Derivation of the NLIE 

Let us recall some basic facts about the NLIE. Consider a 1+1 Minkowski space-time 
discretized along its light-cone directions, with (diagonal) lattice spacing a. The space- 
time has the geometry of a cylinder, of circumference L = Na, where is the number of 
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lattice sites in the spatial direction. On this lattice we consider an inhomogeneous 6- vertex 
model, with an anisotropy parameter 7 and periodic boundary conditions. The Boltzmann 
weights are realized in terms of the quantum R-matrix, in fundamental representation for 
both the horizontal and vertical spaces, of the Uq{sl{2)) afiine algebra, q = —e~^'^. The 
inhomogeneities take the values 6„ = (— 1)"9/2 on the different spatial rows of the lattice, 
so that the unitary time evolution operator can be written as 



U 



-laH 



T 



{2N) ( © 



{©n} 



-(2Af) 



{©n} 



where T^^\i!}\{Qn}) denotes the usual 6- vertex transfer matrix on a square lattice of 
sites (see e.g. [|] for details). 

The diagonalization of this evolution operator is achieved by Algebraic Bethe Ansatz 
with the following set of Bethe equations 



/ 7 
sinh — 

vr 




+ 9 + 


iir 
~2 


7 

sinh — 

TT 




-6 + 


Y 




sinh — 
\ vr 




+ 6- 


in 
~2 


sinh — 

TT 




-e- 


in 
~2 _ 


I 



N 

M sinh - [dj -'&k + ^tt] 

= -n (3-1) 

smh - [^j - - tn] 

TT 

M < is the number of roots in the given Bethe configuration. Notice that the usual 
requirement to have an even number of sites of a Bethe system (in order to be able to 
define an antiferromagnetic vacuum) regards the number 2A^ appearing in the transfer 
matrix t^'^^\ i.e. to the number of sites of the equivalent XXZ chain, and not to the 
number N itself, that can be chosen even or odd without restriction. In general for a 
state with M roots the third component of the spin of the chain is S* = — M. The 
antiferromagnetic vacuum corresponds to the choice of the maximal number of (all real) 
roots (a sort of Dirac sea), i.e. M = N. The energy E and momentum P of a state can 
be obtained from the transfer matrix eigenvalues by use of the formula 



M 



sinh 



7 



Ja(E±P)/2 



IT 



^ 2 



i=i sinh 
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TT 



ITT 



1 + 0(1). 



(3.2) 



The Bethe equations (or rather their logarithms) can be reformulated as quantization 
conditions for the ^j's in terms of the so called counting function, defined as 



M 



Zr,{^) = N[<f)y,{^ + 0) + 01/2(^9 - ©)] - J] 0i(^ - ^k) 



(3.3) 



where 



4>u{^) = i\og 



sinh ■^{inu + i}) 
sinh ^{inu — i)) 
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The branch of the logarithm in 4>u(yi^) is fixed by requiring analyticity and oddness of the 
function in a strip around the real axis. Due to the periodicity of Zn{{}) evident from (|3.3| ), 
the fundamental region to consider is 

27 

which is wider than the analyticity strip. This is the reason for the introduction of the 
so-called second determination of the counting function (cf. [0), to which we come later. 
The roots of the Bethe equations must satisfy the condition 

ZN{^j) = 2nIj , J. gZ + ^, (3.4) 

where 6 = M mod 2 so it takes the values and 1. We return to a discussion of the value 
of S later. 

The physical renormalized continuum QFT on the cylinder is defined by sending a ^ 
and N ^ 00 while keeping L = Na fixed. The renormalized mass scale M. = 4a~^e~® is 
generated by sending the inhomogeneity 6, that can be interpreted as a momentum cutoff 
on the lattice, to infinity together with N, i.e. 

n 1 4^ 
= log — 

(where / is the dimensionless scale I = M.L). In |||] it is shown that there exists a certain 
fermionic operator on the lattice which satisfies equations of motion that on the continuum 
become the Thirring ones. This supports the claim that the continuum limit of this lattice 
theory is describing something related to the sG/mTh models at coupling 

(3 = A/8(7r — 7) 01 p = 1 . 

7 

It is interesting to see to what extent this approach can reproduce the spectrum of sG/mTh 
states. 

The number of roots of the system, which is of order for the low level (non ther- 
modynamical) excitations we are interested in, becomes infinite in the continuum limit. 
The Bethe equations also become infinite in number and are longer useful to describe the 
states. Fortunately, the function Zn{^) has a finite limit, and it satisfies a non-linear in- 
tegral equation (NLIE in the following) already in the lattice, which was constructed some 
years ago |10| directly from the Bethe equations. However, all the derivations reported 
so far, are plagued by one or another imprecision. Therefore we reproduce here a new 
derivation that clarifies some delicate points and also introduces naturally the corrections 
needed for analyticity that were called special roots/holes in ref . [T^ . 



Before attacking this problem, we recall some facts about the solutions of Bethe equa- 
tions that are needed in the following. Real analyticity of the function Zj^{'d) means that 
Bethe roots appear either as real roots or in complex conjugate pairs (with the exception 
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of the selfconjugate wide roots with Qm{} = tt{p+ 1)/2 ). The vacuum state, as said above, 
consists in M = roots, all real. This number becomes infinite on the continuum and can 
be seen as a sort of sea of real roots describing the quantum vacuum of the theory. For 
excited states there can be some points hj on the real axis oft?, where e*^^*^'*^-' = (— l)*^"*"^, 
but that are not roots of the original Bethe equations. These points are known as holes, 
and play the role of particle excitations. Their number Nh is 0(1) and not 0{N) as for the 
real roots, so it remains finite in the continuum limit. The same happens to the number of 
complex roots. It is convenient to distinguish these latter into two classes: the dose roots 
Cj, j = 1, ...,Mc having 

|3mcj| < 7rmin(l,p) , 

and the wide roots Wj, j = 1, M\y with 

• n ^ / ic- I ^ P+ 1 
7rmm(l,pj < |^mtyj| < tt — - — . 

Note that the set of wide roots include, in our conventions, also the selfconjugate roots 
appearing at 



'^mwj = 71- 



p+l 



2 ' 

whose complex conjugate is mapped in the root itself due to in{p + 1) periodicity of the 
counting function. 

2'jv(t9) is real and finite on the real axis and so for every choice for J^'s between the 
minimum and maximum values of Zn{{}) ( p.4| ) has at least on solution, which can be either 
a real root or a hole. Among the real solutions x we distinguish the normal ones with 
Z'^{x) > and the special ones with Z'^{x) < 0. If j/ is a special object (i.e. a special 
root or hole), then its quantum number / is "degenerate", in the sense that there is at least 
another real solution x with the same quantum number: 

ZN{y) = 2nl = Zn{x) , y^x. 

This is implied by the property that the Zisf{{}) function (for N not too small) is "globally 
increasing". Prom the same property it is obvious that the number of solutions correspond- 
ing to a given value / is always odd. 

3.1 Non linear integral equation 

We write the counting function in the following way: 



Nf-. 



Zn{^) = N [01/2(^9 + 6) + 01/2(^9 - 6)] + 5^0i(^ - h 



k=l 



Mc+Mw Mr+Nh 
k=l k=l 



(3.5) 
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where Nh is the total number of holes (normal and special) , Mr the total number of real 
roots (normal and special) and denotes the set of the positions of complex roots {cj , Wj}. 
First we look for an equation for the derivative of the counting function (|3.5|) . Let x 
be a real solution of the Bethe equation. We assume Z'^{x) ^ and put, in a complex 
neighbourhood of x, \ = x + e. Then the expansion 

1 + (-)^e'^^("+^) ^ 1 + (_)5e*^ivW+i.z;^W ^ -itZ'^{x) (3.6) 
has as obvious consequence the identity: 

1 1 (_l)'5e*^iv(^+^)zZ' ix + e) 

^ +... = 7^^(x + e) + ... (3.7) 



A - X e 1 + (-i)'5e*^iv{^+f) 

(the dots are regular terms in A — x). In the following we use the shorthand notation 

> 1 + (_i)<5e^^^W ■ 

One could be tempted to say (as in [l^) that TZ&i^d) is just the derivative of 

log (1 + (-l)V^-W) , 

but this is true only while the logarithm remains in its fundamental domain which is in 
general not true in the lower half plane. This is why we prefer to use the function TZs{'d) 
without expressing it as a logarithmic derivative that could lead to abuse, for example, of 
integrations by parts.0 

The basic trick of the NLIE derivation ||1^, |12], |T3| is to use the Cauchy theorem and 
identity ( p.7|) , to express an analytic function /(x) as 

f{x)=l ^P^={ ^mnsiX), (3.8) 



27ii X — X fr 2ni ' 

where r^, is a curve encircling x anti-clockwise, while avoiding other singularities of the 
denominator i.e. other Bethe solutions (real or complex). This is always possible because 
the solutions are finite in number as long as we are on the lattice with N < oo. Applying the 
same argument applies to each one of the real roots of ( |3.1J ) the last sum in the derivative 
of (|3.5|) can be expressed as: 

Mr+N„ 



E <(^-a^;.)= /^<(^-A)7^,(A). 
k=l -^^ 



This is exactly the origin of unconsistencies in the derivation of [Q. One should either fix the 
logarithm to be in the fundamental domain and pay attention to avoid branch cuts all along the integration 
contour (as we do in what follows) or introduce a "contour logarithm" defined by integrating the function 
TZs along a suitable path taken in the complex plane, but then one needs to pay attention to logarithmic 
identities and the use of integration by parts. 
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X 



X 



X 



— X X a- 



e — X- 



X 



X 



Figure 1: Contour for the integration. The crosses are roots while the circles are holes. 



The sum on the contours was modified to a single curve T, depicted in figure 1, encircling 
all the real roots {x^}, and avoiding the complex Bethe solutions. 

We take T inside the strip — < S3m^9 < 77+ (figure 0), where ri± are positive real 
numbers satisfying 

7T 

< r]+,r]_ < -min{l, p, \'^m^k\ , k = 1, ...,Mc + Mw] ■ 

Without loss of generality we can assume that r]^ = ri_ = 1], and deform T to the boundaries 
of the whole strip. The regions at ±00 do no contribute because 0' vanishes, and so the 
integral can be performed on the lines X = p ±ir], where p is real: 

+00 

—00 
+00 

— CO 

After some algebraic manipulations similar to those in [p!^] one obtains a NLIE for the 
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Z' (i!}) function 



z'm =N 



1 1 



Mc Mw 

- ^ 2nG{^ -Ck)-Yl - ^fc)^^ (3.9) 



k=l k=l 



+253m J dpG{'d — p — i7])TZs{p + 17]) 



where 



+ 00 

G{^) = — [ dke 
2vr J 



TT{p—l)k 



iki) sinh 



2 sinh ^ cosh ^ 



and Gji denotes the so-called second determination of the function G. Using the prescrip- 
tion given in the second determination of any function f{^) is given by 

f^9^ ^ + «7rsign($5m^9)) , p>l, , . 

nf^J/J \ Wsign(S5m^9)) , p<l. ' ^ ' 

Equation (|3.9| ) holds for both positive and negative values of Z'^. We already noted 
that the function IZs is the derivative of 

log [1 + (_i)5e'^^("+''')] , 

but some care is required because of the multivaluedness of the logarithmic function. The 
argument of the logarithm, calculated at a real solution of is (see (|3.6|) ) 



1 + (_i)'5e^^^(-+^+i^) ^ {-ie + 7])Z'^{x) . (3.11) 

It intersects the negative real line which is the cut of the logarithm, if Z'^{x) < 0, which 
corresponds to the case of special roots and holes. If there are no special objects, the 
argument of the logarithm always has a positive real part, without crossing the logarithmic 
cut, thus we can write the integral as: 



+ 00 

— oo 

where log^^j denotes the fundamental determination of the logarithm and the integral runs 
along the real axis. 
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f(y+v) 



f(y) 



f(y-v) 



Figure 2: Values of the function /(x); y is a special object and is a positive real number. 



For real negative w, log(w + iO+) — \og{w — iQ-] 
branch. This is important for special objects yu as evident from (|3.11|) 



2m for any choice of the logarithmic 

Let us denote 



/(x) = 1 + (^—\Ye^^^^^^^'^'> (with X real) and consider the simplest case where there is 
only one special object y {Z'^{y) < 0, Z^iy) = 2ttI and y G M). Observe that \ogp^f{x) 
is discontinuous in y, as shown in figure |^. However, the factor 7ls{x + irj) in ( |3.9| ) is 
continuous everywhere and so it must be the derivative of a continuous function, namely 
of lo 

Spofi^) ~ 2m9{x — y) (plus an arbitrary constant). 
In the presence of Ns special objects we obtain 



Tis{x + irf) 



Ns 



.]_^^5g^2Jv(a;+^^7) 



) - ^ 2'Ki6{x 



Vk) 



k=l 



In the following, we omit the label FD and log will always denote the logarithm in its fun- 
damental branch. It is very important to note that the special objects are not independent 
degrees of freedom unlike the holes and complex solutions that are fixed "a priori": the 
special objects appear when the derivative of at a root or at a hole becomes negative. 
Now we can come back to ([3^ ) and write it with the explicit inclusion of special objects 



N 

Ns 



1 



+ 



Nh 



k=l 



cosh(?9 + 6) cosh(?9 - 6) _ 

Mc 



(3.12) 



k=l 



k=l 



k=l 



J dp ^ 



The odd primitive of G^-d) 
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& oo . , 7r{p-l)k 

//• JL smn 

d,G{x) = y (3.13) 

n _„ 2 sinn cosh — 

oo 2 2 

has the curious property that it equals the logarithm of the soliton-soliton sG scattering 
amplitude Sj^^^-d) as remarked in ||l^]. This is much more than a coincidence and as we 
shall see in section || that it is central in reproducing the sine-Gordon scattering amplitudes. 

Integrating in ^9, we obtain the fundamental non linear integral equation for the counting 
function: 



Z^{^) = 2N arctan — ^ + ^ ;^(^ 2 ^ - Vk)- 

k=l k=l 

Mc Mw 

-^X{'& - Ck) -^X{'& -Wk)n+ (3.14) 

k=l k=l 

oo 

+2$5m j dpG{d-p-i7]) log (1 + (-1)V^'^(''+^'')) + C . 



The integration of the convolution term is possible because G vanishes exponentially at 
infinity. C is the integration constant and must be determined by a comparison with the 
asymptotic values of Zn calculated by the definition (p.5|). 



3.2 Continuum limit 

To perform the continuum limit we send A^, B — > oo in the following way: 



B ^ loe 



ML 



as explained in H]. To obtain the energy and momentum eigenvalues in the continuum 
limit requires a continuum counting function that we define as 



Z{^) = lim Zj^i'd) . 

N oo 



We obtain by performing the limit on ( p.l4| ) 



Z{^) = /sinhT? + gi'dl^j) + 2Qm j dxG{'d - x - ir]) log (l + (-1)^6^^^^+^'')) + C , 

— oo 

(3.15) 
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where I = A4L is the dimensionless volume parameter. The function g{'d\'dj) is the so- 
called source term, composed of the contributions from the holes, special objects and 
complex roots which we call sources and denote their positions by the general symbol 
{^j} = {^k , Vk , Ck , Wk}and takes the form 

Nh Ns Mc Mw 

k=l k=l k=l k=l 

The integration constant gets redefined by the asymptotic value of source terms which 
disappear at the two infinities in the continuum limit. The value of the constant C can be 
obtained using the fact the function Z as defined by ( p.5| ) is asymptotically odd modulo 
27r. The integral term goes to if ^9 — ±oo, while the source terms are odd and so we 
derive that C* = modulo 27r and so we are free to choose (7 = 0. The only problematic 
point arises when there are selfconjugate roots, for which we will work around the problem 
by slightly changing the definition of Xii with respect to the one given in fl^. From the 
formula (|3.10| ) in the repulsive regime we obtain that the source term for a selfconjugate 



root (at position w with ^mw = 7r(j9+ l)/2) is 




cosh^(2(?9-5Rew) + 
"cosh^(2(^?-3fJew) - 




In choosing the logarithmic branches, we require that this function be continuous, but then 
it cannot be odd. Therefore we modify the source term for this case and define it to be 



cosh^(2('i9-3?ew)+i7r) 
^ cosh ^ (2(^? - 3fJe w) - m) 



This is equivalent to shifting Z by an odd multiple of tt. The nice feature is that we can 
then set C = for all source configurations. This redefinition will come up again in section 
|5|, where it will allow us to write the result for the attractive and the repulsive regime in a 
unified form. The notation xii will be used without change but with the above redefinition 
understood where necessary. 

The equation reported in has the form 



Z{d) 



oo 

lsiYi\i^ + g{^\^j)-i j dxGid-x-iri) log (l + + 

— oo 

oo 

+i J dxG - X + IT]) log {{-lY + e-'^^""-'"^) . 



(3.16) 



This equation differs from ours by the position of the term (—1)'' in the second logarithm. 
This difference, however small it appears, has profound implications. First of all, contrary 
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to the claims made by the authors of ||T3|, for 5 = 1, the second term of ( p.l6| ) does 
have logarithmic branch ambiguities even in the absence of special objects; it is enough 
to examine its behaviour substituting Z{i!}) = Zsinh^? which is the limit of the counting 
function for very large I. 

The second point is that with our form of the equation the convolution term decays 
exponentially fast with I. This is important because in the large I limit the source terms 
allow us to reconstruct an asymptotic Bethe Ansatz based on the S-matrix of the continuum 
theory (see section || for details). With the equation (|3.16|) the contribution coming from the 
convolution term decays as a negative power of / if 5 = 1. However, from general principles 
we know that the corrections to the asymptotic Bethe Ansatz come from correlation effects 
which should decay exponentially with the distance in a massive field theory. 

Finally, taking the rj ^ limit one can reconstruct the original Bethe equations from 
( p.l4| ). For that reason one has to examine the function 



QnW = lim - log — — T — .y . . . 



With our formulas we obtain 



Qj^(^) = (Zjv(^9) + nS) mod 27r 



(3.17) 



which has jumps of 27i exactly where the value of {} coincides with the position of a real 
root or hole. By a partial integration one obtains (5-function terms and gets back to the 
original definition of the counting function (p.5| ). 

On the contrary, using the formula ( p.l6| ) one arrives at 



which is wrong for the case 6 



X X + (—lYe'^^^^'^^'^'^^ 

lim - log - — —J .„ . . 

ri->o I i—iy + e-*^Jv(i?-«??) 



Z]^(J}) mod 27r 



1 because the jumps are not at the positions of real 



roots/holes. The authors of ||T^ use the formula (|3.17|) which is the correct one (cf. eqn. 
(4.14) in their paper), but it does not follow from their form of the NLIE. 

The continuum limits for energy and momentum are the similar to the ones exposed in 
|[T3||. The expression for energy, subtracting the bulk contribution is: 



E — E] 



bulk 



Nh Ns 

= M ^ cosh hj — 2M ^ cosh yj 

Mc Mw 

-Ai cosh Cj + M. ^^(cosh Wj)// 

oo 

-M j ^2^m[sinh(x + ir/)log(l + (-l)V^("+'''))] 



(3.18) 
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while for the momentum we get 



P = M ^ sinh hj - 2M ^ sinh yj 

Mc Mw 

—M. sinh Cj + A4 ^^^(sinh Wj)n (3.19) 



oo 



-M j ^2^m[cosh(x + zr7)log(l + (-l)V^("+*''))] , 

— oo 

where the values of hj, Cj, Uj and Wj are fixed by the quantization conditions 

The input required to determine a state is given by the number of holes, close and wide 
pairs and their quantum numbers Ij. It is not necessary to specify the special sources: since 
Z'i^) is dominated by the term Zcosh^9 for large /, it is always positive and so the special 
sources disappear. As a consequence, to specify a state one can give only the normal 
sources at far infrared. The need for the special sources can be detected by gradually 
decreasing the value of / and looking at the sign of the derivative of the counting function 
at the positions of the real roots/holes. 

We remark that there exists a useful relation between the number of different types of 
sources, which in the continuum looks like [|l^ 



Nh - 2Ns = 2S + Mc + 2e{Ti - 2-i)Mw , (3.20) 

which we call the counting equation. The analogous relation on the lattice differs by a 
contribution from sources that disappear in the continuum limit. One can define 

NH,eff =Nh- 2Ns , 

which is in fact the total number of the physical solitonic particles (i.e. solitons and 
antisolitons) in the given state and is independent of /, which we proceed to explain. 
Suppose that J is a Bethe quantum number corresponding to some special objects and let 
us denote by {zj} the set of real numbers for which Z (zj) = 2ttI. First note that since Z is 
globally increasing the number of these positions is always odd. For large enough / the set 
{zj} has only one element Zi, which is normal but can change its nature when decreasing 
/ together with the appearance of new positions {zj : j ^ 1}. All the new positions which 
are generated in this way are holes and satisfy the following properties (for illustration see 
Figure 

1. If for a given / Zi is normal then half of the other {zj : j 7^ 1} are special, while the 
other half is normal. For that reason the part they contribute to Nh is the double of 
their contribution to Ns and so they do not change Nn^eff- Note that special holes 
have to be counted both as special objects and as holes (and therefore contribute to 
the source term with two pieces). 
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e 



Figure 3: The typical behaviour of the counting function Z when a) z\ is normal and when b) 
z\ is special. 

2. If z\ is special, then the number Ns of the special objects in {zj : j ^ 1} is related 
to the number Nh of holes in the same set by Nh = ^{Ns + 1). Together with the 
assumption that Zi is special we obtain a vanishing contribution to NH,eff again. 



4 Scaling functions in the infrared limit 

In this section we analyze the scaling functions coming from the NLIE in the IR limit. 
For / oo the convolution term in the equation becomes small (of order O (e~^)) and 
can be dropped. The surviving part can be seen as a dressed higher level Bethe Ansatz 
giving constraints on the asymptotic states. In particular, the scattering physics can be 
reproduced in this limit and the S-matrix can be read out of the NLIE. This analysis was 



performed in ||13| for pure hole states; here we extend it in order to extract results on the 



quantisation rules of complex roots which will be crucial in the sequel. 



4.1 States with holes only 

As our first example, let us take states with only holes in the sea of real Bethe Ansatz 
roots. In the IR limit we can write 

Z{d) = lsmh^ + J2x{^ - hj) , Z{hj) = 27rlj. (4.1) 

i=i 

Since the function x can be written as 



x{^) = -t\ogSX+{^), (4.2) 

where S^^{'&) is the soliton-soliton scattering amplitude in sine-Gordon theory, the equa- 
tions (|4.1|) describe the (approximate) dressed Bethe Ansatz equations for states containing 



19 



only solitons (or antisolitons). The IR limit puts no restriction on the quantum numbers Ik 
of the holes; this seems to be a very general feature that remains true even when complex 
roots are allowed. 

4.2 Neutral two-particle states 

Let us now extend our investigation to situations with complex roots and consider the 
two particle states in more detail. Forgetting for the moment the breathers, we have 
to consider the two-soliton states. The soliton-antisoliton come in doublets so there are 
four different polarizations for two particle states of which the ss and ss have topological 
charge Q = ±2. These two states are expected to have exactly the same scaling function for 
energy and momentum because the sG/mTh theory is charge conjugation invariant. There 
are two different situations for the neutral ss state, which have spatially symmetric and 
antisymmetric wavefunctions (denoted by (ss)+ and (ss)_, respectively). To separate the 
symmetric and antisymmetric part one simply has to diagonalize the 4x4 SG two particle 
S-matrix and see that it has 2 coinciding eigenvalues (equal to e'^^^^\ corresponding to ss 
and ss), and two different eigenvalues in the Q = channel. 

Now we proceed to demonstrate that the IR limit restricts the possible quantum num- 
bers of the complex roots. To simplify matters we consider only the repulsive regime 
p > 1. 

In the repulsive regime p > 1, the scaling function (ss)_ is realized as the solution to 
the NLIE with two holes (at positions 'i?! ^) and a complex pair at the position p±ia . In 
the IR limit we have 

Zi-d) = /sinh^9 + x(^-'i9i) +x(^-^2) -x(^-p-«cr) -x(^-p + ^cr), 
Z{p ± ia) = 2nl^ . 

The quantization condition for the complex roots explicitly reads (we write down the 
equation only for the upper member of the complex pair, since the other one is similar) 



I sinh(p + ia) + x{p + i(T - i^i) + x{p + icr - ^2) - x(2«(x) = 27r/+ (4.4) 

Now observe that as Z 00, the first term on the left hand side acquires a large imagi- 
nary part, but the right hand side is strictly real. The imaginary contribution should be 
cancelled by some other term. The function x is bounded everywhere except for isolated 
logarithmic singularities on the imaginary axis. For the cancellation of the imaginary part 
the argument 2ia of the last term on the left hand side should approach one of these singu- 
larities (similarly to the analysis in TEA 0). In the repulsive regime, taking into account 

that for a close pair a < tt, the only possible choice for a is to approach ^is / 00. The 

soliton-soliton scattering amplitude has a simple zero at = ivr with a derivative which we 
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denote by C (the exact value does not matter.) To leading order in I, the cancellation of 
the imaginary part reads 



I cosh p + 3f?e log C (^0- - ^ j = , 



(4.5) 



from which we deduce 



TT 



a — 



exp (— / cosh (p)) , 



(4.6) 



so the imaginary part of the complex pair approaches its infrared limit exponentially fast. 
This approach is modified by taking into account the finite imaginary contributions coming 
from the source terms of the holes and from the convolution term. These contributions 
lead to corrections of the order e~' and so they modify only the value of the constant C. 
For the real part we get, again to the leading order 



^ex {p + ^\- ^i) + 3^ex (p + ^\- ^2) = 27r/+ 
It can be shown that (in the repulsive regime p > 1) 



(4.7) 



sinh - (i •§ 

2 



e (^) ^ Uex {& + il) = -\ log ^' j 



sinh - (i— -\- •&] 
p \ 2 / 



(4.8) 



where the branch of the logarithm is specified by ^ (0) = and continuity. ^ is an odd 
monotonous function bounded by 



m\ < ie(oo)i 



^{p - 1) 

2p 



(4.9) 



from below and above. This means that for any allowed value of the real position of 
the complex pair is determined uniquely and that 



p — 1 



2p 



(4.10) 



and since in the repulsive regime p > 1, the only possible choice is — 0. Then the 
solution for p is 



(4.11) 
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Figure 4: The analytic structure of the counting function if o" > ^. The logarithmic cuts are 
indicated with the wiggly lines. The one lying in the upper half plane originates from the root in 
the lower half plane and vice versa. We also indicated the value of the discontinuity across the 
cuts. 

so it approaches the central position between the two holes (the corrections to this asymp- 

totics are also exponentially small for large /). In fact, for the symmetric hole configuration 

Ii = —I2 we expect 'di = —'62 and p = to be valid even for finite values /. 

vr 

However, the above derivation is valid only if cr < — and so we do not cross the boundary 

vr 

of the analyticity strip of the x(^) function which is at 3m ^ = it. If a > — we have to 



use the second determination, which leads us to a different conclusion. Equation (|4.4| ) can 
be written in the form 

/ sinh(p + ia) + x{p + - di) + x{p + icr - 1^2) - x(2^cr)// + x(2«o- - in) = 2nlc , 

where 



= X(^) + - ^tt) = -i log 



sinh ( 

p 



sinh 



The conclusion that a approaches | exponentially fast as Z 00 is unchanged, but now 
the real part of x(2«o") is not zero but instead ivr (modulo 27r). We choose the branches of 
the logarithm in such a way that 3?e x(±2i(T) = ivr (the structure of the cuts is displayed 
in figure ^). In this way we obtain that 

and the asymptotic value of p is just as before. 

The S-matrix of the (ss)_ configuration can be obtained by substituting the asymptotic 
values 



P 



a 



71 

i— 
2 



(4.12) 
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into the expression for Z{'d) and considering now the quantisation rules of the holes. We 
obtain 



Z{^i) = lsmh{A) -Hog S. {^1-^2) = 27r/i , (4.13) 
and a similar equation for the second hole, with 



S4^) = h^Sttm^ (4-14) 



sinh 
sinh 



2p 



which is the correct scattering amplitude for the antisymmetric configuration of the soliton- 
antisoliton system. The amplitude has poles at {} = i7i{l — 2kp) , k = 1,2, . . . , correspond- 
ing to breathers with mass 



m2k = 2M. sin kpn 



The equation ([4.1 3|) describes an approximation to the full NLIE valid for large I which 



is called the dressed or asymptotic Bethe Ansatz. There are three levels of approximations 
to the scaling functions: the full NLIE (which is in fact exact), the higher level Bethe Ansatz 
(HLBA) obtained from the NLIE by dropping the convolution term and the asymptotic BA. 
The difference between the BA and the HLBA is that while the former keeps the complex 
pair in its asymptotic position, the latter takes into account the corrections coming from the 
fact that a changes with / (to first order as given in ([1.6| ) ). However, since the convolution 
term is of the same order in / as the dependence of a derived from the HLBA, the HLBA 
is not a self-consistent approximation scheme. Therefore we are left with the exact NLIE, 
which is valid for all scales and with the BA as its IR asymptotic form. 

The above conclusions about the (ss)_ state an be extended into the attractive regime 
as long as p > ^. At the point p = \, however, the pair situated at the asymptotic 
position ±z| hits the boundary of the fundamental analyticity strip situated at inp. This 
phenomenon has a physical origin: this point is the threshold for the second breather bound 
state which is the first pole entering the physical strip in the (ss)_ channel. To continue 
our state further, it requires to go to configurations having an array of complex roots of 
the first type (cf. ||13|)- Without going into further details, let us mention only that any 
such array contains a close pair plus some wide pairs depending on the value of p. We 



return to comment on this later in subsection [5.3.1 . 

The symmetric state (ss)+ can be obtained from a configuration with two holes and 
one selfconjugate wide root. Considerations similar to the above lead to the scattering 
amplitude 
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cosh 



2p 



cosh 



2p 



stm 



(4.15) 



which agrees with the prediction from the exact S-matrix. The poles are now at 
27r(l — {2k + l)p) , k = 0,1, . . . and correspond to breathers with mass 



fn2k+i = 2M. sin 



{2k + l)p7r 



This configuration extends down to p = 1, where the selfconjugate root coHides with the 
boundary of the fundamental analyticity strip at in. This is exactly the threshold for the 
first breather which is the lowest bound state in the channel. For p < 1 the 

state contains a degenerate array of the first kind (cf. [0]), which consists of a close pair, 
a selfconjugate root and some wide pairs, again depending on the value of p. 

We remark that it is easy to extend the above considerations to a state with four holes 
and a complex pair. In this case the only essential modification to the above conclusions 
is that the limit for the complex pair quantum number becomes 



< 



p 



IT 



when a < —, and 



< 



p 



p — 1 



(4.16) 



p 



(4.17) 



TT 



when ^ since now there are four x sources coming from the holes. In the range 
1 < p < 2 this gives the same constraints 



0, 



a < 



TT 

2 

TT 



=F- , a > — 
^2 ' 2 

as for the (ss)_ state (for p > 2 it allows for more solutions, which however we will not 
need in the sequel). 



5 UV limit and kink approximation 

5.1 Calculation of the ultraviolet conformal weights 

We now examine the ultraviolet limit of the states described by the NLIE in order to 
compare it with the known facts about the UV limit of the sG/mTh theory, outlined in 
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subsection |0. We follow closely the approach described in detail in However, our 



NLIE differs from the one proposed in ||T3l as we explained in section ||. Furthermore, we 



do not fix (5 to be S" mod 2 as the authors of [O] do but instead leave it as a free parameter 



taking the values and 1. For these reasons, it is worthwhile to review the main points of 
the derivation. For simplicity we first suppose that there are no selfconjugate wide roots, 
for which we had to redefine the source term in the repulsive regime (see section |^). 

From conformal perturbation it is known that the leading ultraviolet behaviour of 



energy and momentum are given by (c.f. (6.S)): 



E{L)= -— (c-12(A+ + A_)) + ... , 

27?^ (5.1) 
P{L)= ^(A+-A_) + .... 

Therefore, for the UV comparison we have to determine the contributions to the energy 
and momentum which behave as 1/L when L 0. Consider e.g. the contribution of a 
hole positioned at h to the energy which is mcosh/i. To get a 1/L contribution in the UV 
limit we have to require 

2 

h = finite ± log y , (5.2) 

which means that we are principally interested in the holes whose position has a logarith- 
mically diverging term as I ^ 0. Similar argument can be applied to complex roots. 

The behaviour of the sources for Z — can be classified into three possibilities: their 

position can remain finite (we call them "central"), or they can move towards the two 
2 

infinities as ±logy ("left/right movers"). For sources in all three classes we can introduce 
the finite parts 'df'^ of their positions -dj by extracting the divergent contribution: 



{^,} ^ Uf±\ogj,^'A . (5.3) 



2 



We denote the number of right/left moving resp. central holes by and similarly we 
introduce the numbers A^^'", M^'° and M^'°. 

In the UV limit it is expected that the NLIE can be split into three separate equations 

..e .... ... .e ..... >o. l To ..e 

observe that Z has an implicit dependence from / which we make manifest writing Z = 
Zi^-d, I). Let us define the so-called "kink functions" in the following way: 

Z±{^) =\imZ (^^±\og'^, . (5.4) 

The term "kink" has its origin in the fact that the asymptotic form of the function Z{i!}) has 
two plateaus stretching between the central region and the regions of the left/right movers. 
The functions Z±{{}) describe the interpolation between the plateaus and the asymptotic 
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behaviour of Z{{}) at the corresponding infinity. Note that if there are no central objects 
then the two plateaus merge into a single one stretching from the left movers' region to 
the right movers' one. 

We introduce the quantities S^, S^, by the definition: 

= l[iV±'° - 2Np' - Mp' - 2Mp'9{n - 2j)] 



in analogy with the counting equation ( p.20|) for S, except for the fact that 5* is always 



integer, whereas can be half-integer; moreover, in some cases, can be negative. They 
are interpreted as the spin of the right/left moving and central solutions. It is obvious that 

s = S+ + S~ + 5°. 

We also denote 

7r/2-7 



Xoo = X[ + 00) = 71- 



TT — 7 



For convenience we switch to 7 as the coupling parameter in our UV calculations. Substi- 
tuting the expressions ( |5.3| ) and (|5.4| ) into the NLIE (p.l5| ) and keeping only the leading 
terms as / — >^ we obtain for the left/right movers: 

00 

Z±{^) = ±e^^ + g±{^) + j dxG{^-x)Q±ix) , (5.5) 

—00 

where the source term g±{i!^) is given by 

g±{{}) = ±2xoo {S - S^) + 27r/± +Y,x{^- hf) - 

iV± M± M± 



(5.6) 



i=i j=i i=i 



where is an integer depending on the configuration of the wide roots (and also on the 
regime of coupling taken), whose explicit expression is not important in the sequel. The 
quantization conditions read 

i^f) = 2nlf . 

The equation for the central region can be obtained by taking / ^ in (|3.15|) : 



00 

Zo{^) = goi^) + j dxG{i9 - x)Qo{x) 
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together with 



" H 



j=i j=i j=i 



Zo(^°) =27r/°, 

where Zq{X) = hm^^o Z{X, I). 

The functions Q±fi{x) are connected with the corresponding Z±fi{x) by the usual ex- 
pression ( p.lTp : 

Q.(x) = - log ^ ^ ' -^<QAx)<^, 

where o" = ±, and rj takes an infinitesimal positive real value. 

In the following we will often need the asymptotic values of and Q^j. From ( |5.5| ) we 
obtain for Z±: 

Z+(+oo) = +00 = -Z_(-oo), (5.7) 
the corresponding asymptotic value for Q± is: 

Q±(±oo) = 0. (5.8) 
For the other limiting values of Z± we obtain from (|5.5| ) the plateau equation: 



Z±{toc) = g±{Too) + ^Q^iToo) , (5.9) 

TT 

while from the definition of Q we have 

Z^{^) = Q^(^) + nS + 27ik . 

From (|5.(j| ) we obtain 

g±{Too) =±2{S- 2S^) Xoo + 2iTk^ , 

where again are some integers depending on wide roots. This leads to 

Q^i^oo) = ±2(7r - 27) {S - 2S^) - 2(7r --f)i6 + 2k±) , (5.10) 

k± are integers determined by the condition — vr < Q+(— 00) < tt. Its presence allows us 
to absorb the integers k^. Including selfconjugate wide roots can take 
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half-integer values. By a long but elementary analysis one can see that either both of them 
are integer or half-integer depending on whether the total number of selfconjugate roots is 
even or odd, respectively. This means that in the context of the UV analysis we can define 
an effective value of 6 

Seff = 5 + mod 2) = 5 + (number of self — conjugate roots mod 2) . 

In the repulsive regime there is no ambiguity in the choice of k because 4(7r— 7) > 27r, which 
is larger than the width of the fundamental strip of the logarithm. In some cases compatible 
with the counting equation apparently there is no solution for the plateau equation ( |5.9| ). 
This problem can be overcome by observing that some special sources (roots/holes) are 
required as we will demonstrate on explicit examples later. 

Considering now the region 7 > 7i/2, there can be more than one different possibilities 
to choose k. In ||T3| it is shown that this phenomenon is related to the fact that wide roots 
become excitations independent of the holes, which is manifestly evident from the fact that 
they no longer contribute to the spin of the state S and that the asymptotic value of their 
source contribution xii in the attractive regime is modulo 27r. 

For the central part we obtain the following equations: 

Zo(±oo) = ±2(5 - 2S^)xoo + 2nk^ + ^Qo(±oo) , 

TT 

that coincide respectively with equations for 00) and Q_(+oo). This implies that 

Qo(+oo) = Q+(-oo) , Qo(-oo) = Q-(+oo) , 

corresponding to the formation of the plateau regions as described above. 

Now we are ready to perform the ultraviolet limit on the expressions of the energy and 
momentum. To this end we substitute the ultraviolet behaviour of the sources into the 
energy and momentum expressions and retain only terms of order j^. The kinetic terms 
give contributions that look like 

^e^^^. (5.11) 

The integral terms must be calculated separately. Their form in the energy and momentum 
expressions is, respectively: 

00 00 

-M f ^smhxQix) = -M f ^J[e"Q(x) -e-^'Q(x)], 
J 2t: J 2tt 2 

—00 —00 
00 00 

-M [ ^coshxQix) = -M [ ^l[e^Q{x)+e-^Q{x)], 
J 27r J 2n 2 

—00 —00 

where the contribution depending on is called kink contribution, and the other term 
is the kink contribution. For the computation observe that, in the limit of very small /, 
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it is possible to write an expression for Q: 

Q{x, I) r^Q-(^x + log +Q+(^x~ log - Qo{-x) + q{x, I) , (5.12) 

where the function q{x, I) vanishes in the / limit. 

Substituting (|5.12|) in the integral form of the "+" kink contribution we obtain 



oo oo 

— oo — oo 

where the symbol = means that only the terms in 1/L are retained. Similarly we get for 
the "-" kink term 

oo oo 
— OO — oo 

Now we are ready to express energy and momentum in a way dependent only on quantities 
which are finite in the UV limit. For the conformal energy (scaling weight) we obtain 



(T=± \ j = l j = l j = l j = l 



ZTT 



and the conformal spin is given by 



^ - ^ = lElE - 2 E -^'"^ - E + E " (^°"' ) 

\j=i 



II 

tT=± \ j = l j = l j = l j = l 

OO 



27T 

(T= 
OO 

dx 



27r 



Using the method explained in [1^, it is possible to give a closed expression for the sums 
and the integrals. The computation is simple, and we perform it explicitly only for the 
case of holes. The quantisation condition for a hole hf , using the kink equation (|5.5|) and 
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the oddness of the x function allows us to write: 

k=l 

















j=i \ 


k=l 




K 




+ 








k=l 






'^H poo 





(5.15) 



The contribution from the other sources can be evaluated in a similar way. Summing 
together all the expressions obtained in this way the terms containing x cancel completely 
because of the oddness of the function x- In the integral term we make the substitution 



27iG{x) = x{x) . (5.16) 



The result is: 



In 



c 

~ 24 


oo 


1 


f dx 




J 2^ 



where 



/oo 
dxGi^ ~ x)Q±ix) , (5.17) 
-oo 



w ' 



(5.18) 



and we denote 



is an integer or half-integer depending on the configuration of the wide roots, which is 



best calculated case by case as it has already been noted in 1 13 



Using the "kink lemma" proved in ||I3|, the integral term can be evaluated: 



dx 
2^ 



V±{x)Q±{x) 



± 



12 87r(l-7/7r) 
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Our final expression for the conformal weights is 

Let us discuss certain properties of the formula ( |5.19|) . Using the relation 



1-2 1 



TT 2i?2 

one can prove with simple algebra that the conformal weights have the form 



where 



m = 2S 



and 



For a consistent interpretation of the weights A± ( |5.20|) should have the form (p?^). Match- 
ing the terms we get the condition 

nl = n^_ . (5.22) 

The second condition is that one must be able to choose the sign of n in such a way that the 
left and right descendent numbers (obtained by combining N± with a term depending on 
n and m) are integers. Since the contribution comes exclusively from the last term in 
( ^.19| ), the one-plateau systems satisfy trivially the condition (|5.22|) , but even in that case 



it is not clear whether the descendent numbers are integers. At the time of this writing we 
cannot exclude the possibility that there exist some states for which the conformal weights 
cannot be interpreted within the framework of c = 1 CFT. In the numerous examples we 
have calculated so far we have not found any such source configuration. If it happens for 
some configuration of the sources then the corresponding scaling function must be excluded 
from the spectrum. 



From the second line of (5.21) we obtain a useful relation 



2n± = (5 + Msc) mod 2 i.e. 2n± = Sefj mod 2 , (5.23) 

where we used that 25 and 45^ are even. 

We also drive attention to the fact that we have no proof that the scaling functions 
obtained by the method of NLIE span the complete space of states with even topological 
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charge. This is an extremely difficult problem due to the following two circumstances: (1) 
the dependence of the UV conformal weights from the parameters of the source configu- 
ration is rather complicated and (2) to obtain the allowed values of the complex roots one 
has to examine the IR limit as in section ^. It is also clear that the same state can be 
realized by different root configurations depending on the value of the coupling constant 
(see e.g. the {ss)± states in subsection ^]2|). The investigation of this question is out of 
the scope of this paper. 

Using the relation (|5.23| ) it can be determined what quantisation rule one has to impose 



on the real roots/holes if one aims to get solutions that in the UV limit reproduce the states 
of the sG/mTh theory. The rule is 

Seff = , 

which is different from the rule 6 = S mod 2 advocated by Destri and de Vega in ||T3|. Their 
result is based on the reasoning that in the continuum limit the length of the spin chain 

goes to infinity, there is no loss generality in taking even. In fact, this assumption 
fixes the rule 6 = S mod 2. In our approach, however, we do not fix N to be even on the 
lattice, which allows us to leave the parameter 6 free and to fix it later when we examine 
the UV limit of the states described by the NLIE. For states without selfconjugate roots 
our rule simply means that the real roots/holes are always quantised by half-integers. 

An interesting phenomenon is that the conformal weights obtained depend only on very 
generic features of the source configuration such as the asymptotics of the left and right 
moving sources and the total spin. This means if we have a certain source configuration, 
we can add new sources separately to the right and the left moving part in such a way 
that they are separately neutral (i.e. do not change and S~). In this way we do not 
change the primary weight, but generally we increase the term — 21 f — — and 
so create descendents of the state we started with. An example: states which have S = 
and = are all descendants of the vacuum, however complicated their actual source 
configurations are. 



In subsection |3]^ it was claimed that the best way to specify a state in the NLIE 
framework is to give the corresponding source configuration for large / since then the 
special sources are absent. In the UV asymptotic calculations, however, we work just at 
the opposite end of the range of scales. Therefore a different method to obtain the special 
sources is necessary. The approach we apply is to find a regime in the coupling 7 where 
the plateau equation (|5.9|) admits a solution without special sources and then determine 
the solution in other regimes by analytically continuing the values of Z-|-(=i=oo) in 7. This 
will allow us to determine the special sources and will provide a way to fix the ambiguity of 
the plateau solution (|5.10|) in the attractive regime as demonstrated on explicit examples 



in the following. 

Finally, let us remark that putting 6 = S mod 2 our results for the conformal weights 
are identical to the ones obtained by Destri and de Vega (formulas (8.21), (8.22) and 
(8.23) in ||l3)- We now go on to give a detailed examination of states with simple source 
configurations in order to clarify the consequences of the general result ( |5.19| ). Some of 
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the examples correspond to configurations whose scahng functions will be compared to the 
TCS data in section ^ 



5.2 Examples of states without complex roots 
5.2.1 The vacuum state 

Starting from a lattice with 2N sites one finds that the antiferromagnetic ground state 
which has spin 5* = 0, when written in terms of Bethe vectors depends on M = roots of 
the Bethe equations, all of which are real. We expect that this ground state corresponds 
to the vacuum of the field theory. However, there are two possibilities: one for N even and 
the other for N odd, corresponding to 5 = or 5 = 1. As we show in the sequel, only one 
of these states can be identified with the vacuum for a local field theory having a c = 1 UV 
limiting CFT. The UV dimensions of the vacuum are A''^ = 0. Choosing 5 = we obtain 

c=l, 

as in pop . 

For the choice 5 = 1 case we obtain an interesting result. From (|5.10|) we derive 



Q+(-oo) = -2(7r-7)(l + 2fc), (5.24) 
which admits solution only in the attractive regime when 

Q+(-oo) = ±2(7r -7) . (5.25) 
For the plateau value of Z we obtain 

Z+(-oo) = ±(7r-27). 

The value is not determined uniquely due to the contribution of wide roots, but it can be 
fixed using information from the repulsive regime. 

In the repulsive regime, the real root at the origin is a special one with Z'iGi) < 0. The 
sources are composed of two normal holes and the special root at the origin, all of which 
has Bethe quantum number zero. One of the holes moves to the left and the other one to 
the right and so = ^, while the special root is central. This allows for a unique solution 
of the plateau equations: 

Q+(-oo) = 27 , Z+(-oo) = -TT + 27 . 

Let us suppose that the counting function Z is analytic as a function of 7 which is plausible 
because examining the NLIE for this state we find that all terms are analytic in the 
coupling. This determines which branch to choose for the plateau values in the attractive 
regime. We obtain the final result (for both the attractive and the repulsive regime) 

A± - _ 2^ - ^ - A± 
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a) 



b) 



Figure 5: The UV behaviour of Z for a) 7 > ^ and b) 7 < |. In b), the root at the origin turned 
into a special one and two normal holes appeared at the points where Z crosses the horizontal 
axis with a positive derivative. 

which are the conformal weights of the vertex operators V(±i/2,o)- The actual UV limit 
can be a linear combination of these operators; in any case, it is not contained in the UV 
spectrum of the sG/mTh theory. 

Let us analyse the solution briefly to see the mechanism responsible for the appearance 
of the special sources in the repulsive regime. Since the source term is zero in the attractive 
regime, we expect that the Z function is odd in with the behaviour depicted in figure |^ 
a). Passing through the free-fermion point 7 = | to the repulsive regime, we observe that 
the plateau values change sign and the derivative of the function at the origin becomes 
negative (part b) of figure |^ so the root located there turns into a special one. At the 
same time, however, there appear two new holes, one left and the other right moving, with 
Bethe quantum number 0. This phenomenon can be thought of as the splitting of the root 
at the origin into a special root and two holes and is one of the general mechanisms in 
which the special sources are generated. Another mechanism will be explored when we 
examine the UV limit of the soliton-antisoliton states. 



5.2.2 States containing two solitons 

As we have explained, for a generic configuration there are two possibilities, corresponding 
to the half number of sites being even and odd. For the case of two solitons (i.e. Nn^eff = 
Nh — ^Ns = 2), the simplest possibilities are: 

1. Two holes quantised with J12 = ±|, 6 = 0. Using the results of the previous 
subsection yields: 



corresponding to the operator V(o,2), which is the UV limit of the lowest-lying two- 
soliton state as it can also be seen from the TCS data. If instead of the minimal 



4(l-7/7r) 



1 




0,2 ' 
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choice Ji 2 we take a nonminimal one Ji = /+, /2 = — /- with I± G Z+ — ^ we obtain 



2' 

1 

2 ' 



which corresponds to descendents of V(o,2)- This is a general phenomenon: the "min- 
imal" choice of quantum numbers yields the primary state, while the nonminimal 
choices give rise to descendents. 



2. Two holes quantized with J12 = ±1, 5 = 1, as proposed in ||I3|. In the repulsive 
regime a special root is required, as in the case of the 5 = 1 vacuum state. The result 
(in both regimes) is: 

^" = i + irz^ = ^"A2 = Ar/2, + i- 

This state is a linear combination of a_iV(i/2,2) and a-iV(_i/2,2) which means that it 
is not contained in the local operator algebras of sG/mTh theories. 

3. Two holes quantized with /i = 0, /2 = ±1, 5 = 1. Consider in detail the case 
of /2 = +1 since the other one is similar. Suppose that the hole with /i = is 
a left mover and the other one is a right mover (the other possibilities lead to a 
contradiction). We find a solution to the plateau equation only in the attractive 
regime: 

Q+(-oo) = ±2(7r-7), Q_(+oo) = ±2(7r-7). 

In the repulsive regime the hole with quantum number becomes a special hole y and 
emits other two ordinary holes each quantised with zero. We obtain for the plateau 
values 

Q+(-oo) = Q_(+oo) =27. 
The conformal weights turn out to be: 

1 (7A)' 1 / 1 

41 - 7/7r 2 \2R 

These are the conformal weights of the vertex operator V(i/2,2)- Performing a similar 
calculation for I2 = —1 we obtain the weights of V(-i/2,2)- 
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5.2.3 Generic number of solitons 



As the calculation proceeds in complete analogy with the two-soliton case, we only give 
a summary of the results. The conformal families which we obtain depend on how many 
holes move to the left and to the right. The primaries are obtained for the minimal choice 
of the hole quantum numbers; increasing the quantum numbers we obtain secondary states, 
as it was pointed out in the case with two holes. The states we obtain are in the conformal 
family of a vertex operator \4,m, with 

2 



m = NH,eff , n G Z + - 



As a consequence, all states with 6 = are contained in the UV spectrum of the sG/mTh 
theory while the ones with 6 = 1 are not, in agreement with the relation (|5.23|) . The 



complete expression for n is somewhat complicated, but we give it in the case of symmetric 
( = Njj = Nu,effl'^ ) configurations: 

n = , 5 = 0, 
n = ±-, 6 = \. 



5.3 Examples with complex roots 
5.3.1 Two holes and a close pair 

We now turn ourselves to the state with two holes and a close pair (in the repulsive regime) , 
which describes the antisymmetric soliton-antisoliton two-particle state. We recall that 
according to the results from the IR asymptotics, the quantum numbers of the complex 
roots can be or ±|. Therefore, we are left two possibilities: either the holes are quantised 
with integers or with half-integers. 

Let us start with the half-integer choice 5 = and suppose that one of the holes is right 
moving (which is the case if its quantum number l\ = /+ is positive) and the other one is 
left moving (i.e. I2 = I- < 0). Using the general formula (|5.10| ) and 5*^ = |, 5 = 0, by a 
simple calculation we obtain 

Q+(-oo) =Z+(-oo) = -2(7r-27), 

which is valid for f < 7 < f (the complex pair remains central) . The other plateau values 
follow by oddity of the function Z . For the other half of the repulsive regime we have 
to include two new normal holes, one moving to the left and the other to the right, and 
two special holes that remain central (a justification for this will be given shortly) so that 
5± = 1. We then find 

Q+(-oo) = 47 , Z+(-oo) = -2(7r - 27) , 
which is valid for 7 < |. In both cases the result for the UV weights turns out to be 

A^ = l- -±/±-i = ± /± - - = Af±i 0) ± /± - - . 

TT ^ 2 2i?2 ^2 ^^^'"^ ^ 2 
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The primary state is obtained for the minimal choice /+ = — /_ = | and coincides with 
a linear combination of the vertex operators V(±i,o) which is correct for the state to be 
included in the spectrum of the sG/mTh theory and agrees with the behaviour of the 
{ss)_ state observed from TCS. If both of the holes move to the right or to the left, we 
obtain descendents of the identity operator, i.e. states in the vacuum module of the UV 
CFT. 

Figure ^ shows how the special holes are generated in this case. Starting from 7 > f 
and decreasing the value of 7 the UV asymptotic form of the counting function varies 
analytically. Since the real roots/holes are quantised by half-integers, they are at positions 
where the function Z crosses a value of an odd multiple of vr. As the plots demonstrate, 
the behaviour of Z is in fact continuous at the boundary 7 = ^: it is our interpretation 
in terms of the sources that changes, exactly because we try to keep the logarithm in the 
integral term of the NLIE in its fundamental branch. The price we pay is that we have to 
introduce two new normal holes (one moves left and the other one moves right) and two 
special holes which are central. 

Let us comment on the extension to the attractive regime. The complex root config- 
uration changes only by the possible presence of wide pairs (see subsection |4^ ), which 



however have no effect on the plateau values ( |5.10| ) as their contribution can be absorbed 
in redefinition of the values k±. The other effect the wide roots have is to shift the terms 
S-t by multiples of 27r, but this affects only the integers N± in ( p.2q ). Therefore we again 
obtain states which are descendents of V(-|-i 0). 

For the state with integer quantisation of the holes, we only give the result. If one of 
the holes is moving left and the other one to the right, we obtain 

1 / 7\ ^ 1 



i.e. some linear combinations of descendents of the vertex operators V(±i/2,o)- If both of 
the holes move to the left or to the right, we obtain other descendants of the same primary 
states. Again, this excludes the integer quantised states from the spectrum of sG/mTh 
theory. 
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5.3.2 Two holes and a selfconjugate complex root 

Taking the integer quantised state, the plateau values are given by the same formulas 
as for the case with the close pair above, since the plateau equation is identical with 
6eff = 6 + 1 mod 2 = 0. The only difference is that the numbers g^j^g in ( |5.1^ ) take a 
nonzero value 



If the one of the holes is a right mover, while the other one is a left mover, the conformal 
weights turn out to be 

A± = 1 - ^ ± /± - 1 = A(±i,o) ± /± - 1. 

TT 

If the two holes move in the same direction, we again obtain secondaries of the vacuum 
state. 

Concerning the extension to the attractive regime, one can again check that the plateau 
solution remains unchanged for the root configuration described in subsection therefore 
the conformal family to which the state belongs remains the same, similarly to the case of 
the state. 

One can ask what happens if we quantize with half-integers. The result is, similarly 
to the case with the complex pair, that we obtain descendants of the operators V(-|-i/2,o)- 
Therefore we conclude that in this case the integer quantised configuration must be ac- 
cepted, while the half-integer one is ruled out, in full accordance with the rule 5e// = 0. 



5.3.3 Four hole examples 

We do not discuss all the possible examples of four holes with complex roots. Instead, we 
concentrate only on two cases that will come up in the comparison we make with TCS 
data in section |^. 

The first example is four holes with a complex pair. We choose the half-integer quantised 
prescription and the hole quantum numbers Ji = —I2 = | and I3 = —1^ = |, which is the 
minimal choice with two left moving and two right moving holes (the nonminimal choices 
lead again to descendents). Performing the calculation we obtain that 

A± = ^- ^ + 2 = — + — - + 5.26 

47r(7r-7) 2 2R^ ^ ' 

The operators having these conformal weights are the following: 

a-2V(2,i) ,a_ia_iV(2,i) ,a-2V(2-i) and a_ia_iV(2 -1) . 

The UV limit of the scaling functions corresponds to some linear combination of these 
states. This is also confirmed by comparison with TCS data. The integer quantisation 
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leads again to states not present in the UV spectrum. Considering four holes with a 
selfconjugate root in the middle, we arrive at similar result, with the difference that this 
time the integer quantisation leads to the correct weights ( ^.26| ), while the half-integer 
prescription yields states which are not part of sG/mTh Hilbert space. 

6 Truncated Conformal Space at c = 1 

Before presenting the numerical results, let us give a brief description of Truncated Con- 
formal Space (TCS) method for c = 1 theories. The TCS method was originally created 
to describe perturbations of Virasoro minimal models in finite spatial volume 0. Here we 
present an extension of the method to study perturbations of a c = 1 compactified boson, 
more closely the perturbation corresponding to sine-Gordon theory. 

The Hilbert space of the c = 1 theory can be split into sectors labelled by the values of 
P and Q, which are quantised by integers. In the numerical computations we shall restrict 
ourselves to the P = sector (it is not expected that any relevant new information would 
come from considering P ^ 0). The TCS method consists of retaining only those states in 
such a sector for which the eigenvalue of He ft is less than a certain upper value E^ut, so 
the truncated space is defined as 



nTcs{s,m,E,ut) = {1"^) : P |^) = s |^) , Q |$) = m |^) , i/cFT |^) < ^c«t |^)} • (6.1) 

For a given value of s, m and Ecut this space is always finite dimensional. In this space, 
the Hamiltonian is represented by a finite size matrix whose eigenvalues can be computed 
using a numerical diagonalization method. The explicit form of this matrix is the following: 



H= — {Lo + Lo--I + X^—tzkB 1 , (6.2) 



where Lq and Lq are diagonal matrices with their diagonal elements being the left and 
right conformal weights, I is the identity matrix, 



is the scaling dimension of the perturbing potential V and the matrix elements of B are 

B^,^ = i ($1 V(i,o)(l, 1) + V^(o,i)(l, 1) 1^) . (6.4) 
We choose our units in terms of the soliton mass A4 which is related to the coupling 



constant A by the mass gap formula obtained from TEA in 19|: 
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(6.5) 



where 



2T{h/2) 

rv 



/ _ \ 

TT 



2-h 



7rr(l - h/2) 



2r 

V \2-hJ \A-2h) J 



h 



(6.6) 



In what follows we normalize the energy scale by taking Ai = 1 and denote the dimen- 
sionless volume AiL by I. For numerical computations, we shall use the dimensionless 
Hamiltonian 



H _2n ^ c - '2-/. 



h=T:^ = — {^o + Lo-—I + B I . (6.7) 



The usefulness of the TCS method lies in the fact that it provides a nonperturbative 
method of numerically obtaining the spectrum (the mass gap, the mass ratios and the 
scattering amplitudes) of the theory. Therefore it can serve as a tool to check the exact 
results obtained for integrable field theories and get a picture of the physical behaviour even 
for the nonintegrable case. The systematic error introduced by the truncation procedure 
is called the truncation error, which increases with the volume L and can be made smaller 
by increasing the truncation level (at the price of increasing the size of the matrices, which 
is bound from above by machine memory and computation time). 

Let us make some general remarks on how the TCS method applies to c = 1 theories. 
First note that the Hilbert space (even after specifying the sector by the eigenvalues of 
P and Q) consists of infinitely many Verma modules labelled by the quantum number n. 
At any finite value of Ecut only finitely many of such Verma modules contribute, but their 
number increases with E^ut- As a result one has to deal with many more states than in the 
case of minimal models. The results of TCS are supposed to approach the exact results in 
the limit E^ut oo, and the convergence can be very slow, while the number of states rises 
faster than exponentially with the truncation level. The perturbing operator has scaling 
dimension which ranges between and 2, becoming more relevant in the attractive regime, 
while the number of states corresponding to a given value of E'eut becomes larger as moving 
towards p = 0, which affects the convergence just the other way around. 

Generally, the energy of any state goes with the volume L as 

k=l 
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where (A^i) are the left (right) conformal dimensions of the state in the ultraviolet 
limit, B is the universal bulk energy constant (the vacuum energy density) and the infinite 
sum represents the perturbative contributions from the potential V. 



The bulk energy constant has also been predicted from TBA and reads [IS 



B = -]tan(^) (6.9) 



4 V 2 

(the same result was obtained from the NLIE approach in ||lOl). This is a highly nonanalytic 
function of p and it becomes infinite at the points where p is an odd integer. In fact, at 
these points there is a value of k for which k{2 — h) = 1, and (^) oo. The infinite 
parts of B and Ck i'^) exactly cancel, leaving a logarithmic (proportional to / log/) and a 
finite linear contribution to the energy, by a sort of a resonance mechanism. All of these 
"logarithmic points" are in the repulsive regime. However, due to UV problems in the 
repulsive regime we are not able to check numerically the logarithmic corrections to the 
bulk energy. 

The origin of UV divergences can be understood from conformal perturbation theory 
(CPT). It is known that when the scaling dimension h of the perturbing potential exceeds 1, 
CPT suffers from ultraviolet divergences which should be removed by some renormalization 
procedure. The TCS method is something very similar to CPT: it operates in the basis 
of the UV wave functions as well, but computes the energy levels using the variational 
approach and therefore could be called "conformal variation theory" (CVT). As a result, 
we expect that there could be UV divergences for the range of couplings where h > 1 which 



is exactly the repulsive regime p > 1 ||20|. The numerical analysis has in fact shown that 
in the repulsive regime the TCS energy eigenvalues did not converge at all when increasing 
the truncation level. 

Fortunately, there exists a way out: since we expect to find a sensible quantum field 
theory when the UV cutoff is removed, it should be the case that the relative energy levels 
£m{L) = E^,{L) — Ejjac{L) converge to some limit. This is exactly the behaviour that we 
observed. Consequently, in the repulsive regime one can only trust the relative scaling 
functions produced by the TCS method, while in the attractive regime even the absolute 
energy values agree very well with the predictions of the NLIE [|l^ (including the predicted 
bulk energy constant ( |6.9|) , which is completely analytic for p < 1 and thus logarithmic 
corrections are absent as well). 

Our numerical results show that the smaller the value of p is the faster the convergence 
is (with the understanding that in the repulsive regime by convergence of TCS we mean 
the convergence of the energies relative to the vacuum). On the other hand, even in the 
attractive regime the convergence is so slow that to get reliable results (which means errors 
of order 10~^ — 10~^ for the volume / ranging from to somewhere between 5 and 10) one 
has to work with dimensions around 4000. This means that the TCS for c = 1 theories is 
far less convergent than the one for minimal models (in the original Lee- Yang example the 
authors of took a 17 dimensional Hilbert space (!) and arrived to very accurate results). 
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7 Comparison of the NLIE to TCS 



In this section we give the numerical comparison between the NLIE solutions and the 
TCS data. In the IR limit of the NLIE gives back the scattering theory of the sG/mTh 
model. In this limit, the scaling functions are described by the asymptotic BA, which for 
two-particle states can be written as follows 

I sinh A - i log S{A - ^2) = 27r/i, 

lsmh^2-i^ogS{^2-A) = 27Tl2, (7.1) 

where ^1^2 are the rapidities of the two particles and the energy of the state relative to the 
vacuum is given by 

^hhi^) = "^1 coshi?! + m2 cosh^92, (7.2) 

with the understanding that one has to substitute for t^i 2 the solutions of the equations 
( fni) . The corrections from the NLIE to the IR asymptotic Bethe Ansatz of ( fTll) decrease 
exponentially with /. This is also expected on physical basis: the theory has a finite mass 
gap and hence these finite size corrections should decay exponentially with the volume. 
Therefore to see a substantial effect of the NLIE as compared to the asymptotic BA one 
has to go down to small values of the volume. This is also advisable for the comparison 
with TCS, because the truncation effects are smaller if one is closer to the UV region. 

On the other hand, the UV dimensions computed by the "kink approximation" were 
shown to be correct. Therefore the comparison between TCS and NLIE serves in fact for 
verifying that the scaling functions obtained from the NLIE interpolate correctly between 
the UV and IR limits in the intermediate range of scales. Since the accuracy of TCS is 
getting worse for larger volumes and, on the other hand, the iterative solution of the NLIE 
becomes problematic for small volumes, it is at intermediate scales that a useful comparison 
can be made. 

To have a comparison between TCS and the NLIE, we have to get numerical predictions 
for the scaling functions from the NLIE. Before presenting the results we briefiy discuss 
the method used for the numerical solution of NLIE. 

We recall that the NLIE (|3.15|) contains a source term g{'d\'dj) which is specified by 



giving the root/hole structure of the given state and depends on the positions of the 
holes, special roots/holes and complex roots. These positions in turn are fixed by the 
Bethe quantisation conditions. The coupled nonlinear equations coming from the integral 
equation supplemented with the quantisation conditions can be solved numerically by an 
iterative procedure. First one chooses a starting position for the positions of the sources 
'dj, which we set to the asymptotic position for the complex roots and to zero for all hole 
positions. Then one iterates the integral equation (using fast Fourier transform to evaluate 
the convolution) to update the counting function Z. Using this new Z, an improved 
determination of position of the -dj can be obtained, and is fed back into the integral 
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equation for a new iteration cycle. The process is repeated until the solution is found to 
a prescribed precision ( usually 10~^ ). For large /, the source term dominates, while the 
correction coming from the integral term is exponentially small. Therefore we expect that 
the further we go to the IR regime the faster the iteration converges which does in fact 
hold. Hence it is preferable to start iterating at the largest desired value of / and decrease 
the volume gradually, always taking as a starting point at the next value of the volume 
the solution found at the previous value. 

We do not present results for pure hole states which have already been discussed in a 
previous paper 0], both in the attractive and the repulsive regime. Instead we concentrate 
on the new results obtained for states containing complex roots and in the repulsive regime. 

Let us start with the {ss)± states in table [l| and ||, respectively. The NLIE results can 
be considered to be exact to the accuracy showed, therefore the total deviation must be 
attributed to truncation effects. Hence we expect that it shows a growing tendency with 
the volume /, which is in fact the case (for the lowest lying {ss)± states the NLIE and 
TCS predictions intersect each other around I = 7.5 and / = 6.0, respectively, but this is 
a numerical artifact due to truncation errors from TCS). The order of the deviations is 
consistent with what we expect from measuring the accuracy of the mass gap from TCS. 
It also gets larger if we consider a state of higher energy which is expected since we get 
closer to the truncation level. 

It is also possible to consider four hole states with complex roots (see table |). In this 
case we consider the lowest example of four holes with a selfconjugate wide root resp. with 
a close pair, with the prescription of integer resp. half-integer quantisation. 

The reason to confine ourselves to the lowest choices for the quantum numbers is related 
to the problem of identifying the states, since the higher we go in energy the denser the 
spectrum gets. Given that the results of TCS are approximate, it becomes very hard to 
identify the points in the TCS plot which form a line corresponding to a given energy 
level. This job can be done by using BA combined with other information. For example, 
it can be observed that the states (ss)+ and (ss)_ always come in pairs close to each 
other, the symmetric state being higher. These energy levels cannot cross each other 
because the sG/mTh theory is parity conserving and these states have opposite parity. 
This information is enough to identify the first two such pairs up to / = 10. The four- 
particle states in table |^ lie in the sector with topological charge Q = 2, where the first 
few energy levels are soliton-soliton two-particle states. Eliminating these states from 
the TCS spectrum (using asymptotic BA or even full NLIE predictions) facilitates the 
identification of the required lines. The reason why the data presented in table |^ only go 
up to Z = 5 is that after this point the TCS points corresponding to the given state cannot 
be identified unambiguously. For the same reason we do not present four hole states with 
higher quantum numbers. Although there are TCS data points which agree with the NLIE 
predictions within the estimated truncation errors, the density of the TCS points is so high 
that it makes impossible to identify the correct ones. 
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h = -/2 = 2 


I 


TCS 


NLIE 


Dev. (%) 


TCS 


NLIE 


Dev. (%) 


0.5 


15.4665 


15.4152 


0.33 


40.3852 


40.2870 


0.24 


1.0 


8.0919 


8.0485 


0.54 


20.3794 


20.2325 


0.73 


1.5 


5.7245 


5.6729 


0.9 


13.7695 


13.5837 


1.4 


2.0 


4.5887 


4.5325 


1.2 


10.5064 


10.2886 


2.1 


2.5 


3.9358 


3.8782 


1.5 


8.5782 


8.3347 


2.9 


3.0 


3.5186 


3.4617 


1.6 


7.3148 


7.0501 


3.8 


3.5 


3.2317 


3.1775 


1.7 


6.4283 


6.1475 


4.6 


4.0 


3.0236 


2.9734 


1.7 


5.7750 


5.4820 


5.3 


4.5 


2.8662 


2.8211 


1.6 


5.2757 


4.9737 


6.1 


5.0 


2.7430 


2.7038 


1.4 


4.8826 


4.5749 


6.7 


5.5 


2.6437 


2.6113 


1.2 


4.5658 


4.2546 


7.3 


6.0 


2.5619 


2.5367 


1.0 


4.3054 


3.9931 


7.8 


6.5 


2.4932 


2.4756 


0.7 


4.0878 


3.7762 


8.3 


7.0 


2.4344 


2.4248 


0.4 


3.9034 


3.5941 


8.6 


7.5 


2.3834 


2.3820 


0.06 


3.7451 


3.4394 


8.9 


8.0 


2.3385 


2.3456 


0.3 


3.6077 


3.3068 


9.1 


8.5 


2.2987 


2.3143 


0.7 


3.4874 


3.1921 


9.3 


9.0 


2.2629 


2.2872 


1.1 


3.3811 


3.0923 


9.3 


9.5 


2.2305 


2.2635 


1.0 


3.2864 


3.0047 


9.4 


10.0 


2.2010 


2.2428 


1.8 


3.2017 


2.9275 


9.4 



Table 1: Numerical comparison of the NLIE prediction and the TCS data for the first two (ss)+ 
states, at p = 1.5 . 



44 





2 


h = -l2 = - 

2 


I 


TCS 


NLIE 


Dev. (%) 


TCS 


NLIE 


Dev. (%) 


0.5 


15.0922 


15.0605 


0.21 


40.2909 


40.1949 


0.24 


1.0 


7.6401 


7.5916 


0.64 


20.2522 


20.1042 


0.73 


1.5 


5.2475 


5.1891 


1.1 


13.6247 


13.4355 


1.4 


2.0 


4.1198 


4.0565 


1.6 


10.3511 


10.1282 


2.2 


2.5 


3.4900 


3.4264 


1.8 


8.4171 


8.1667 


3.1 


3.0 


3.1009 


3.0408 


1.9 


7.1500 


6.8771 


4.0 


3.5 


2.8428 


2.7889 


1.9 


6.2620 


5.9710 


4.9 


4.0 


2.6621 


2.6163 


1.8 


5.6076 


5.3039 


5.7 


4.5 


2.-5298 


2.4934 


1.5 


5.1082 


4.7952 


6.5 


5.0 


2.4292 


2.4032 


1.1 


4.7158 


4.3967 


7.3 


5.5 


2.3502 


2.3352 


0.65 


4.4000 


4.0779 


7.9 


6.0 


2.2864 


2.2827 


0.17 


4.1410 


3.8184 


8.5 


6.5 


2.2337 


2.2415 


0.35 


3.9250 


3.6041 


8.9 


7.0 


2.1892 


2.2085 


0.88 


3.7422 


3.4251 


9.3 


7.5 


2.1509 


2.1818 


1.4 


3.5858 


3.2738 


9.5 


8.0 


2.1175 


2.1599 


2.0 


3.4503 


3.1450 


9.7 


8.5 


2.0879 


2.1417 


2.5 


3.3318 


3.0343 


9.8 


9.0 


2.0614 


2.1264 


3.1 


3.2274 


2.9385 


9.8 


9.5 


2.0374 


2.1134 


3.6 


3.1347 


2.8551 


9.8 


10.0 


2.0156 


2.1023 


4.1 


3.0517 


2.7821 


9.7 



Table 2: Numerical comparison of the NLIE prediction and the TCS data for the first two (ss)_ 
states, at p = 1.5 . 





4 holes (-2,-1,1,2) + s.c.w. 


4 holes j 


'3 1 1 3\ 

, 2' 2'2'2j+^-P- 


/ 


TCS 


NLIE 


Dev. (%) 


TCS 


NLIE 


Dev. (%) 


0.5 


51.2829 


51.2224 


0.12 


50.7801 


50.7101 


0.14 


1.0 


26.2457 


26.1581 


0.33 


25.6191 


25.5111 


0.42 


1.5 


18.0383 


17.9318 


0.59 


17.3879 


17.2521 


0.79 


2.0 


14.0138 


13.8965 


0.83 


13.3812 


13.2257 


1.2 


2.5 


11.6458 


11.5258 


1.0 


11.0490 


10.8823 


1.5 


3.0 


10.0998 


9.9802 


1.2 


9.5416 


9.3709 


1.8 


3.5 


9.0162 


8.9016 


1.3 


8.4957 


8.3272 


2.0 


4.0 


8.2162 


8.1101 


1.3 


7.7318 


7.5706 


2.1 


4.5 


7.5924 


7.5086 


1.1 


7.1518 


7.0015 


2.1 


5.0 


7.1110 


7.0380 


1.0 


6.6972 


6.5609 


2.1 



Table 3: Numerical comparison for the first four hole states containing a selfconjugate wide root 
resp. a closed pair, with the four holes having the quantum numbers indicated in the table. 
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8 Conclusions 



In this paper we have studied how the nonlinear integral equation deduced from the light 
cone lattice model of ||^ describes the excited states of the sine-Gordon/massive Thirring 
theory. We can summarize the new results as follows: 



1. We have presented a new derivation of the fundamental NLIE from the light cone 
lattice which avoids the problems originating from the multivaluedness of the complex 
logarithm function and gives us a correct form of the NLIE. 

2. By examining the infrared limit of our equation we have shown that it leads to the 
correct two-particle S-matrices of sine-Gordon theory. 

3. By computing the UV conformal weights from the NLIE we have shown that it is 
consistent with the UV spectrum of sG/mTh theory only if we choose the parameter 



6 (i.e. the quantisation rule) in a way different from the proposal of ||l3l. The new 



rule is that 6 should be equal to the number of the self-conjugate roots modulo 2. 

4. We have verified the predictions of the NLIE by comparing them to results coming 
from the TCS approach. 



In view of our results we think that the philosophy behind the NLIE has to be changed. 
Namely, the NLIE describes nothing else than a collection of scaling functions including 
ones not realized in the quantum field theory we wish to study. We have seen that the 
spectrum provided by the NLIE is much larger than necessary for sG/mTh, and the correct 
approach is that from the variety of scaling functions predicted by the NLIE one has 
to select the ones corresponding to the Hilbert space of the quantum field theory to be 
described. The key condition in performing the selection is the locality of the resulting 
operator algebra (which we established at the UV fixed point, using the results of ITtI ). 
With this point understood, our results provide a strong evidence that the NLIE does in fact 
describe the finite volume spectrum of sG/mTh theory (for even value of the topological 
charge) . 

Let us discuss briefly the relation between the TBA approach to excited states 
and the NLIE approach. The TBA approach is based on an analytic continuation of the 
vacuum TBA equation to the complex plane of the volume parameter /. By encircling 
certain branch points in that plane, it is possible to obtain the equation describing the 
evolution of some excited state scaling function and repeating the procedure other states 
can be found. At present this is the most powerful approach known to obtain excited state 
scaling functions in perturbed minimal models. The main drawback of the method is that 
it is difficult to compute the equation for a generic excited state. 

The NLIE approach has the advantage that it can give the equation for all of the excited 
states in closed form. The price we pay is that it is a complex nonlinear integral equation 
while the TBA equations are real, therefore its structure is much more complicated due to 
the fact that it encodes a lot of information in compact form. Another limitation is that 
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it requires the knowledge of an integrable lattice regularisation. On the other hand, this 
has a positive consequence too: we do not put in by hand the model to be described as 
in TBA, but it emerges during the derivation (see the comments about the appearance of 
the sine-Gordon phaseshift in section ^ after equation ( p.l3| )). 
Finally, let us point to some open questions: 

1. Is the set of scaling functions provided by the NLIE complete i.e. can we find to 
every sG/mTh state (of even topological charge) a solution of the NLIE describing 
its finite volume behaviour? This can be called a "counting problem". The main 
difficulty is that the structure of the solutions is highly dependent on the value of 
the coupling constant - to see that it is enough to consider e.g. the appearance of 
special sources. There should also be dependence on the coupling related to physical 
reasons (e.g. breather thresholds - cf. subsection |4.2|). 



From the form of the source terms in the NLIE it seems likely that the excited state 
equations can be obtained by an analytic continuation procedure analogous to the 
one used in TBA to obtain the excited state TBA equations. Certain features of 
the arrangement of the complex roots in the attractive regime and their behaviour 
at breather thresholds also point into this direction. This an interesting question to 
investigate because it can shed light on the organization of the space of states and 
can lead closer to solving the counting problem described above. 



3. It was conjectured in ||12| that a so-called a-twisted version (a la Zamolodchikov 
[0) of the NLIE can lead to excited state scaling functions of $(1,3) perturbations of 
minimal models which could provide a link between the NLIE and the TBA approach 



to excited states. In |112| it was found that the vacuum scaling functions obtained in 
this way agreed with the ones predicted by TBA within a very small error. Therefore 
studying the a-twisted NLIE is clearly an interesting problem. 

The NLIE technique so far has been limited to states with even topological charge. A 
work on the description of the odd sector and of the sine-Gordon - massive Thirring 
difference is under way, the results of which will be published soon pT|. 



The second and third points are also important because their investigation may lead closer 
to understanding the relation between the TBA and the NLIE approaches. It is quite likely 
that establishing a connection between the two methods would facilitate the development 
of both and may point to some common underlying structure. 
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